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The  deterministic  prediction  of  both  transient  and  CW  field 
coupling  to  large  complex  electrical  systems  poses  such  a 
formidable  problem  that  other  approaches  to  predicting  the 
coupling  may  be  useful.  In  this  paper,  we  analyze  the  contin- 
uous wave  coupling  to  a complicated  electrical  system  in 
terms  of  random  small  dipole  interactions  in  the  low  frequency 
limit  (wavelengths  >>  system  components'  sizes).  Both  random  . 
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coupling  to  the  incident  wave  and  random  interactions  among 
the  dipoles  are  considered.  The  variables  being  randomized 
are  the  incident  direction  and  polarization,  the  sizes  and 
orientations  of  the  dipoles,  the  mutual  coupling  strengths, 
and  the  lumped  load  impedances.  The  resulting  normalized 
current  distributions  are  shewn  to  be  insensitive  to  the 
details  of  the  model  except  at  the  extremely  low  and  high 
percentiles . _ 

The  magnetic  dipole  case  is  investigated  in  detail.  Its 
resulting  induced  current  distribution  roughly  resembles, 
but  is  not,  a log-normal  distribution  with  a standard 
deviation  in  the  vicinity  of  ^ 6-7  dB.  This  result  provides 
insight  into  some  recent  measurements  obtained  for  EMP 
transient  field  coupling  to  large  systems. 

An  important  implication  of  the  results  is  that  for  a 
variety  of  complicated  systems,  essentially  consisting  of 
many  small  elements  that  the  coupling  is  dominated  by  low 
frequency  magnetic  fields,  the  central  parts  of  the  induced 
current  probability  distributions  are  similar  and  nearly 
log-normal.  However,  conclusions  based  on  the  extrapolation 
of  log-normality  from  measured  values  near  the  median  to  the 
extreme  percentiles  may  be  susceptible  to  sizeable  errors. 

Presented  in  the  appendices  to  indicate  a larger  scope  of 
the  statistical  approach  are  the  basic  results  for  some 
simple  electric  dipole  cases  and  an  investigation  of  the 
difference  in  effects  caused  by  an  elliptically  polarized 
incident  wave  versus  a linearly  polarized  one. 
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SECTION  1 

INTRODUCTION  AND  SUMMARY 


Not  only  does  electromagnetic  (EM)  coupling  to  an 
electronic  system  take  place  through  intentional  penetrations 
such  as  antennas  and  waveguides,  but  in  addition  there  usually 
are  present  many  inadvertent  EM  coupling  paths  through  cables, 
apertures,  grounding  loops,  etc.  While  usually  designed  to 
handle  the  normal  signal  and  noise  background  adequately,  an 
unhardened  system  subiected  to  an  unusual  EM  disturbance  such 
as  a radiated  continuous  wave  (CW)  field  or  the  electromagnetic 
pulse  (EMP)  produced  by  a nuclear  detonation  may  be  caused  to 
malfunction  by  spurious  signals  introduced  through  these 
inadvertent  coupling  paths.  Analyzing  and  predicting  inad- 
vertent coupling  for  the  purpose  of  assessing  and  protecting  a 
complex  system  has  historically  been  and  still  is  a difficult 
and  challenging  task. 

While  in  principle  an  arbitrarily  accurate  analysis  of 
the  EM  coupling  can  be  derived  by  solving  the  Maxwell  equa- 
tions in  the  context  of  a boundary  value  problem  [1],  in 
reality  even  for  a relatively  simple  system,  such  a classical 
deterministic  approach  often  demands  more  effort  and  resources 
than  are  available.  To  keep  the  mathematics  tractable,  judi- 
cious use  of  approximations  and  engineering  judgments  are 
inevitably  required.  Even  so,  the  effort  presently  needed 
to  obtain  approximate  deterministic  predictions  for  EM  coupling 
to  complicated  systems  is  still  substantial  (2).  On  the  other 
hand,  when  the  statistical  properties  of  the  coupling  by 
themselves  play  a dominant  role  in  the  problem  being  addressed, 
a much  less  detailed  probabilistic  approach  may  be  adequate. 

In  view  of  these  observations,  a statistical  theory  like  the 
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the  one  presented  in  this  paper  may  provide  a productive 
new  approach  to  EM  interaction  technology. 

The  problem  of  describing  the  EM  interaction  with  the 
many  interaction  paths  and  modes  of  a system  has  some  similarity 
to  the  problem  of  describing  the  behavior  of  a gas  in  a box. 

One  possible  approach  to  analyzing  the  gas  problem  is  to  work 
out  the  interaction  forces  and  then  apply  Newton's  law  to 
solve  simultaneously  for  the  dynamics  of  all  the  gas  molecules 
— a rigorously  correct  approach  which  has  met  with  very 
limited  success.  An  alternate  approach  is  to  take  advantage 
of  the  fact  that  there  are  a large  number  of  degrees  of  freedom 
in  the  motion  of  the  molecules  and  to  formulate  a statistical 
model  of  the  gas  [3].  While  this  statistical  approach  reveals 
much  less  information  than  the  deterministic  one  and  predicts 
only  the  statistical  properties  of  the  gas,  it  has  nevertheless 
proven  to  be  of  considerable  practical  value.  Similarly,  the 
EM  interaction  with  a complicated  system  consisting  of  a large 
number  of  similar  basic  coupling  elements  can  be  modeled 
statistically  to  yield  a number  of  interesting  results. 
Applications  of  these  results  include  estimating  the  distri- 
bution of  EM  coupling  data  to  be  obtained  for  large  systems, 
and  interpreting  the  experimental  results. 

The  fundamental  approximation  in  the  probabilistic 
approach  to  EM  coupling  is  the  hypothesis  of  a basic  element 
of  interaction  whose  parameters  are  statistically  distributed. 

We  take  this  basic  element  of  interaction  to  be  the  small 
dipole  (although  a more  complex  element  could  be  accommodated 
by  the  theory)  and  the  basic  variables  randomized  are  the 
incident  polarization  and  direction,  the  sizes  and  orientations 
of  the  basic  dipoles,  the  mutual  coupling  strengths,  and  the 
lumped  load  impedances.  The  basic  elements  and  variables  are 
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then  aggregated  in  a probabilistically  rigorous  manner  to 
yield  the  final  results.  In  contrast  to  the  deterministic 
approach,  such  probabilistic  models,  when  applicable,  do  not 
become  proportionally  more  complicated  but  should  yield 
increasingly  accurate  results  as  systems  being  analyzed 
become  more  complicated,  if  this  complication  results  from 
the  repetition  of  the  basic  element.  However,  we  emphasize 
that  the  validity  of  the  probabilistic  approach  to  EM  coupling 
should  be  reviewed  for  each  case  to  which  it  is  applied. 

For  complicated  systems  when  the  coupling  is  dominated 
by  low  frequency  magnetic  fields,  a wide  range  of  basic 
coupling  element  parameters  and  statistical  distributions 
produce  a current  probability  distribution  whose  central 
part  is  nearly  log-normal  with  a standard  deviation  of 
about  6 to  7 dB . However,  conclusions  based  on  the  extrapo- 
lation of  the  log-normality  of  the  current  distribution  from 
the  values  near  the  median  to  the  extreme  percentiles  are 
susceptible  to  substantial  errors.  These  results  may  provide 
some  insight  into  similar  EMP  coupling  data  to  large  systems 
[4]  obtained  for  PREMPT  (Program  for  EMP  Testing)  [5]. 

To  indicate  the  scope  of  the  statistical  coupling  method, 
an  investigation  of  the  basic  results  for  some  simple  electric 
dipole  cases  is  also  presented,  as  is  tne  difference  in  effects 
caused  by  the  elliptically  polarized  incident  waves  and 
linearly  polarized  ones.  In  the  following  text,  Section  2 
presents  the  statistical  physical  models.  Section  3 considers 
the  magnetic  dipole  interaction  case,  and  shows  a comparison 
with  PREMPT  experimental  data.  Section  4 presents  results 
for  specific  electric  dipole  cases,  while  Section  5 investi- 
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gates  the  difference  in  effects  as  caused  by  elliptically 
versus  linearly  polarized  incident  waves. 


The  results  obtained  for  subresonant  CW  excitation 
appear  to  be  sufficiently  interesting  to  warrant  investiga- 
tion of  the  response  of  resonant  electric  structures  to 
transient  electromagnetic  waves. 
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SECTION  2 

STATISTICAL  PHYSICAL  MODEL 


2.1  BASIC  MODELS 

Consider  a plane  EM  wave,  with  a linearly  polarized 

magnetic  field  H and  a frequency  w , incident  upon  a thin 

and  small  circular  loop  [6],  as  shown  in  Figure  1.  To 

obtain  a basic  statistical  model,  we  randomize  both  the 

incident  field  and  the  loop  configuration  to  have  the  random 

variables  (0.  , ) for  the  incident  wave  direction,  ¥ for 

in  in 

the  incident  magnetic  polarization,  (0,<5)  for  the  loop's 
normal  direction  n,  and  R for  the  loop's  radius.  Accordingly, 
the  induced  current  has  an  amplitude  [7] 


I = 


; ~ I 

'loop 


% 


^T-UAqJ?  * n 

/-  lu>L  + R , + R. 

w rad  loss 


(1) 


where  u is  the  permeability  of  the  surrounding  medium  (usually 


its  vacuum  value  u ) , 

o 


A the  loop's  area,  L,  its  inductance, 
o w 


and  Rracj  and  R^oss  resPective^y  its  radiation  and  wire-loss 
resistance.  For  a small  and  thin  loop,  the  induced  current 
random  variable  I [8],  from  (1),  can  be  normalized  such  that 
the  assumed  value  i is  given  by 


i = n—  0 < i < 1 

rm  ~ - 


(2) 


with 


| COS0 

cos  ( $ 


cosi p sin9  . 

in 

- . ) + sinij> 

m 


sinQ  [cosi|i  cos9  . 

in 

sin  - 4>i  ) ] | . 


(3) 
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Figure  1.  Orientation  of  the  Incident  Field 
and  the  Loop  with  Normal  n 


Here,  r<r,0<9<ir,  0 < <J> . < 2tt,  0 < 9.  < rr,  0 < <J>. 

— in  — — — in  — — m — — in 

< 2tt  , and  0 < 4i  < 2tt  are  the  values  taken  respectively  by  the 

random  variables  R,  0,  $ , 0.  , , and  V with  their  corre- 

m in 

sponding  probability  densities  denoted  by  pR(r),  pQ  ^(9,$), 

n.  (9.  <b . ) . and  d...  (iM.  In  terms  of  these,  the  normal- 

^0.  $ . m,  in  r 

in,  m 

ized  induced  current  has  a probability  distribution 


F<1oop)(  , , r p (9,0)  d9d*  p (9in.*inl 

J.  ' in,  in 


r < 

- n 


d0ind(f>in  P'V(^)  ^ PR(r)dr 


and  a probability  density 


p(loop) (i) 


d p(loop) 
= dl  F 


{I  < i} 


Equations  (2)  to  (6)  give  the  procedure  to  compute  the 
probability  distribution  of  the  induced  current  from  those 
of  the  randomized  model  parameters. 

For  the  coupling  mechanism  based  on  short  electric 
dipoles,  as  shown  in  Figure  2,  we  can  model  it  similarly. 
For  this  case,  the  induced  current's  amplitude  is  [9] 


M l2,E 


L?2  • 


ln(fj) 


I 


where  L is  the  length  of  the  wire  with  direction  unit  vector 
L,  M the  position  along  the  wire  at  which  the  current  is 

considered,  D the  diameter  of  the  wire,  and  E the  incident 

% 

electric  field.  Compared  with  the  magnetic  coupling  case, 

A 

the  differences  are  that  E,  L,  and  L replace,  respectively, 
H,  n,  and  R and  that  the  additional  random  variables  M and  D 
are  introduced.  From  Equation  6,  the  normalized  induced 
current  assumes  the  values 


nsx  in  b 

1 2 

b in  x 


0 < i < 1. 


(7! 


Here  n is  given  by  (3),  0 £ s £ 1 is  the  value  assumed  by 
the  random  variable  S = M/(L/2)  with  probability  density 
Pg(s),  and  1 <<  x < b is  the  value  assumed  by  the  dimension- 
less length  X = L/D  with  density  pv(x)  where  b is  some 
largest  length  constant.  Finally,  from  (6)  and  (7),  we  have 
the  normalized  induced  currents'  probability  distribution 


Flwire)(I  < i)  - 


f 


sxjn  < i _b_ 

in  x — in  b 


p0  (9,0)  d9d0  p0_  (0in^in) 

m,  in 


d9ind<J)in  dl^  ps(s)  ds  pX(x)  dx  (8) 


and  the  corresponding  density 


p<wire)(i)  = d.  F(wire){1  < .} 


(9) 
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The  above  simple  coupling  models  form  the  basis  of  our 
statistical  coupling  analysis  of  induced  currents.  Before 
specific  cases  are  examined,  we  shall  make  several  observa- 
tions. First,  the  ratio  of  the  current  in  the  magnetic  case 
to  that  in  the  electric  case  is 


R 


I 'magnetic  I , 4R  l"(s) 

n r 2 

1 electric'  ttL  we 


H 

E ' 


(10) 


This,  for  a plane  EM  wave  with  | E/H | = /p/e , reduces  to 

4R  tn(ji)  2R  Un(§) 

R * 2 = 2~2  • " * (ID 

ttw/uT  L it 

Thus,  for  small  and  thin  loops  and  wires  of  comparable  sizes, 

2 

i.e.,  2R  ^ L,  R 'v  X in  (L/D) /2tt  r)  >>  1,  the  magnetic  coupling 
is  much  more  important  than  the  electric  one.  This  leads  us 
to  place  more  emphasis  on  magnetic  field  coupling. 

Second,  some  very  simple  probability  densities  for 
angular  orientations,  to  be  used  later,  are  listed  in  Table  1 
for  convenience. 

Third,  consider  an  observation  that  is  very  simple  but 
useful  in  interpreting  some  of  the  special  results.  For  a 
random  variable  U = f(Y)|cos0|,  where  0 is  the  angle  between 
a fixed  axis  and  an  orientation  equally  probable  in  all 
spatial  directions,  Y is  positive  with  a probability  density 
PY(y),  and  f (y)  > 0 and  is  monotonic  such  that  f 1 (y)  exists. 
Then  the  probability  density  for  U 
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TABLE  1.  SIMPLE  ANGULAR  ORIENTATION  DISTRIBUTIONS 


Description  p (9,<j>) 

d .'P 


Equally  probable 

in 

all 

directions 

sine 

in  space 

4tt 

Equally  probable 

in 

all 

directions 

5(9  - 9 ) 

on  a conical  e = 

eo 

surface 

Fixed  direction  at  e 

i = 

V * = *0 

5(9  - 9q)  6(<f>  - <i>0) 

P0<u>  - / pY(£_1iy>)  [4  f'1(y)]  f (12) 

u 

implies 

up^(u)  = -PY  (f-1  (u))  f"1  (u)]  . (13) 

Thus,  if  f(y)  is  a monotonically  increasing  (decreasing) 
function  of  y,  the  p^Cu)  is  a monotonically  decreasing 
(increasing)  function  of  u.  One  of  the  simplest  examples 
[10]  is  when  f (y)  = y,  in  which  case  U is  the  projected  length 
onto  a fixed  axis  of  a random  length  whose  spatial  orientation 
is  equally  probable  in  all  directions.  Then,  no  matter  what 
probability  distribution  that  random  length  has,  the  Pul1-1) 
is  a decreasing  function  of  u — the  projected  length  is  more 
probable  to  have  smaller  values  than  larger  values. 
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Fourth,  we  emphasize  that  the  simple  models  considered 
are  for  electrically  small  basic  coupling  elements,  loops  or 
wires,  with  the  interactions  among  these  elements  neglected. 
(Such  interactions  may  be  neglected,  for  example,  when  the 
elements  are  placed  far  apart  so  their  interactions  are  small.) 
However,  these  mutual  interactions  can  be  taken  into  account 
and  moderately  change  the  resulting  current  distribution  in 
a predictable  way,  as  will  be  shown  in  the  next  section  and 
the  special  examples. 

2.2  RANDOMIZED  MUTUAL  COUPLINGS  AND  IMPEDANCE  LOADS 

Consider  first  the  mutual  couplings  by  examining  the 
current  induced  on  the  kth  element  due  to  the  incident  field 
and  the  multiple-scattered  fields  produced  by  currents  at 
all  other  elements  [11] : 

- “]JEo + £ Eji  - Ik°> i 1 + Gk'-  (U> 

j^k 


Here,  the  superscript  "o"  stands  for  the  mutual-coupling- 
neglected  current.  Now,  since  the  ^ E.  is  the  sum  of  a 

j^k  3 

large  number  of  fields  of  fully  randomized  current  elements, 
the  central  limit  theory  applies  [12] 


Og,),  a normal  distribution  with  mean 

deviation  a„,.  However,  interested  only  in  the  relative 


Thus  -*■  G ' eN  (uG  1 
and  standard 


G'  * 

amplitude  of  the  induced  current,  we  can  normalize  out  the 


resulting  factor  (1  + u 


G’ 


from  (14)  and  get 


I = I (o) i 1 + G|  (15) 

where  GeN(0,  a ^ , //I  + , ) . Thus  the  standard 

deviation  o of  G measures  the  strength  of  the  mutual 
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coupling,  with  aG  <<  1 representing  very  weak  and  negligible 
mutual  couplings  and  £ 1 representing  strong  coupling. 

The  value  of  o.  is,  of  course,  determined  by  the  more  detailed 
physics  and  geometry  of  the  basic  interacting  elements  and 
cannot  be  derived  from  statistics  alone.  We  shall  treat  it 
as  a parameter. 

From  (15)  and  making  use  of  the  formula  for  the  product 
of  random  variables,  we  immediately  obtain 


where  G'  = 1 + G e N(l,a_).  Thus,  the  inclusion  of  the  mutual 
couplings  among  the  basic  dipoles  gives  rise  to  an  extra 
factor  to  the  induced  currents  as  shown  in  (15)  and  changes 
the  current  probability  density  from  P ^ (i)  to  Pj(i)  accord- 


ing to  Equation  (16).  To  the  mutual-coupling-neglected  and 
normalized  current  I ^ , which  distributes  itself  in  the 
interval  [0,1],  such  an  inclusion  tends  to  enhance  the 
probability  density  for  smaller  i and  to  give  rise  to  a 
prolonged  density  tail  for  i > 1.  For  the  mean  u and  the 
standard  deviation  o^,  a simple  calculation  using  (16)  shows 


= u 


(o) 


VI 


0(.e-1''<2oG) 


*’/ 


1/ aG 


dx 

/'Tn 


-x2/2 


(17a) 
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oG  « 1 


(o) 


[l  * 0(aG)] 


(17b) 


aG  >:>  1 


(O)  ' V?  • t1  + °(4). 


(17c) 


and 


(- 

W 


(O) 


1 + 


(o) 


1 + o 


fl 


, r1/0G  , -1/  ( 2x2 ) 

a e*1/(2aG)  + 2 C dx  e 
G 


i 


/Tir 


- 1 


(18a) 


aG  <<:  1 ( — 

► \ U 


(o) 


(o) 


°(°1) 


(18b) 


0_  >>  1 7T 

► 2 


1 + 


Oj  (o) 
WT(o) 


- 1 + o(—. 

' o‘ 


:i8c) 


Further,  one  can  easily  show  that  the  ratio  £l  + (a^/u^J^j/ 

[l  + (°  (0)/u  (o))2]  ;''s  monotonicaHy  increasing  with  a G> 

Thus,  this  ratio  lies  between  1 and  tt/2  for  all  o_,  a variation 

G 

not  too  sensitive  to  the  mutual  coupling  strength 

Next,  consider  the  effect  of  a randomized  lumped  load 
impedance  Rs  + iXg  in  series  with  the  dipole  elements.  This 
renders  an  induced  current 
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j 


1 


F. 


(19) 


(s) 


= I 


(o) 


5 I 


fo) 


V( 


1 + 


X 

C 

F 


where  Xq  is  the  dominant  reactance  of  the  dipole  element. 
Lacking  the  detailed  information  on  the  distribution  of  the 
load  impedances,  we  make  the  assumption  that  the  probability 
density  associated  with  an  impedance  is  a decreasing  function 
of  the  value  of  the  impedance.  This  is  certainly  true  above 
some  values  of  impedance,  and  considering  the  prevalence 
of  low  impedance  ground  loops  in  most  systems,  that  value 
should  usually  be  zero.  A simple  try  would  be  to  have  a 
Pl/F ( £ ) = qe  ^ ^ for  £ > 1 and  q > 0,  which  includes  the 
case  of  no  load  impedance  as  the  parameter  q •*  °°  and  the  case 

of  a wide  range  of  virtually  uniform  loads  as  q -*  0 , and 

consequently  results  in  a p (f)  = qe+c*  f ^ e for  0 <_  f <_  1 . 

In  view  of  its  having  a strong  zero  at  f = 0 and  being  a smooth 

and  monotonically  increasing  function  of  f in  [0,1]  but  diffi- 
cult to  manipulate  mathematically,  we  can  now  assume  and 
examine  a qualitatively  similar  but  simpler  Pp(f)  to  gain 
insights  into  the  effects  of  such  a distribution--e . g . , by 
choosing 


0 < f < 1 

p„(f)  = (n  + 1)  fn,  (20) 

f 0 < n 


With  this  p_(f),  we  obtain 
F 


P 

I 


(s) 


(i) 


(n  + 1)  in 


n + 1 
x 


(21) 
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and  from  which 


n + 1 
n + 2 


(22) 


(n  + 2 ) 2 1 

(n  + 1) (n  + 3)  (n  + 1)  (n  + 3) 

(23) 


To  effect  a strong  zero  for  the  p„(f  -*■  0)  of  (20),  n must  be 

r 

much  greater  than  1.  Consequently,  from  (21)  to  (23), 
randomized  loads  have  little  influence  on  the  induced  currents' 
statistics,  except  slightly  enhancing  the  near-zero  portion 
of  the  induced  current  probability  density  by  v(n  + l)/n  and 
reducing  the  value  at  its  end  to  zero.  These  results  suggest 
that  it  may  be  reasonable  to  ignore  the  load  impedance  random- 
ization in  our  special  case  investigations. 
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SECTION  3 

SPECIAL  CASES  FOR  MAGNETIC  LOOPS 


In  this  section  we  investigate  the  random  coupling 
results  for  several  cases  of  the  magnetic  loop  model.  This 
could  represent,  for  example,  the  EMP  coupling  with  a com- 
plicated system  consisting  of  many  ground  loops.  We  calculate 
first  the  1^°^,  then  take  into  account  mutua1  coupling  by  the 
method  outlined  in  Section  2.2. 

3.1  VERTICAL  PLANE  ORIENTATIONS 

Consider  first  a conductor  geometry  where  nearly  all 
of  the  loops  are  in  vertical  planes,  such  as  might  be  the 
case  for  a single  large  bay  of  electrical  equipment.  If  all 
the  loops  are  in  vertical  planes,  i.e.,  their  normals  are 
horizontal,  and  have  their  normal  pointed  equally  probably  in 
all  azimuthal  directions,  then  only  the  horizontal  magnetic 
field  component  can  induce  current  in  the  loops  and  all 
horizontal  field  directions  have  the  same  effect.  Thus, 
without  suffering  any  loss  of  generality,  we  can  consider  a 
H^-  polarized  plane  wave  incident  along  the  -z  direction 
(see  Figure  1)  onto  these  loops. 

Now,  if  we  further  assume  that  the  loops'  sizes  are 

2 

distributed  with  equally  probable  areas  between  a.  = Ttr. 

2 11 
to  a^  = r 2 > then  from  (2!  to  (5)  we  obtain  the  probability 

distribution  of  the  normalized  induced  current  amplitudes 


(loop)  U(o)  ,» 
(A, 2d)  I1  - 


r 2 “ 


/ 

v ) 


r dr  d4> 


(24) 
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re  [rL ,r2] 

4>e  [0,2tt] 
r cost)  <_  i r2 

Here,  the  subscript  (A, 2d)  indicates  the  case  of  equally 
probable  areas  of  loops  in  two-dimensions.  The  probability 
density,  from  (24),  is  (see  Figure  3) 

( A - i2  - /a2  - i2  , 

P(^°P)  (i)  = j- 

I ° ( A , 2d)  tt  ( 1 - a ) | ^ ~y 

(25) 


0 < i < g 
a < i < 1 


where  a = r^/r^  <_  1 is  the  parameter  showing  the  range  of 
spread  of  the  loops'  sizes.  Thus,  the  induced  normalized 
current  I ^ has  a most  probable  value  a,  an  average 


(o)  4 (1  + a + a2) 

p (A, 2d)  3tt  (1  + a) 


(26) 


and  a standard  deviation 


(o)  ( 1 + a2 

a(A,2d)  j 4 


1/2 


4(1  + a + g' 

3 tt  ( 1 + g) 


(27) 


If  we  assume  the  loops  have  equally  probable  radii, 
instead  of  areas,  between  r^  and  r2,  we  then  similarly  obtain 
(see  Figure  3) 
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Pj(0)(i) 


NOTE:  P(!°?p)  (0)>P(1°?P)  (0)  FOR  ALL  a 

I 0 (R,2d)  I ^ (A, 2d) 


Figure  3.  Sketches  of  p^]n?P^  (i)and  p^]n?P^  (i),  the 

PU;(A,2d)  rU;(R,2d) 

Normalized  Incident  Current  Probability  Densities 
for  Equally  Probable  Areas  and  Radii  Cases, 
as  Functions  of  i 
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In 


( loop) 

I (o) (R, 2d) 


(i) 


1 + A - i2 


a + /a 


fl  72 
/a  - 1 


7T  (1  - ct) 


iln 


+ A - i‘ 


0 < i < a 


a <_  i <_  1 

(28) 


where  the  subscript  (R,2d)  indicates  the  case  of  equally 
probable  radii  in  two-dimensions.  Its  average  and  standard 
deviation  are  [13] 


(o)  _ 1 + a 

w (R,  2d)  tt 


(29) 


,(o)  _ j 1 + a + 

MR, 2d)  J 6 


1 + a 


-i  2 


1 1/2 


J I 


(30) 


Notice  the  inequalities  uj^2d)  < u|°;2d)  ahd  oj°;2d)  < 

an  intuitively  plausible  result  because  the  equal-probable-area 
case  gives  less  weight  to  small  sizes  than  does  the  equal- 
probable-radius  case.  Also,  (25)  and  (28)  approach  the  same 

expression  2/^tt/I  - i^ ) when  a -►  1,  as  they  should. 


The  a dependences  of  the  averages  and  the  standard 
deviations,  as  well  as  their  ratios,  are  shown  in  Table  2. 

We  see  that  the  standard  deviations,  at  about  0.3,  are  very 

insensitive  to  the  range  of  the  spread  of  x (whether  a wide 

2 2 
spread  over  0 to  -rrr2  or  a narrow  one  clustered  near  rr r 2 ) and 

the  nature  of  the  spread  (whether  equally  probable  in  area 

or  in  radius) . Also  quite  so  insensitive  are  the  ratios  of 


26 


ABLE  2.  AVERAGES  AND  STANDARD  DEVIATIONS  FOR  MAGNETIC  LOOPS  EQUALLY 
PROBABLY  ORIENTED  IN  VERTICAL  PLANES  WITHOUT  MUTUAL  COUPLING 


. rl 
a = — 

r2 

u(0) 

u(A.2d) 

_ ( o ) 

0 ( A , 2d ) 

0 A. 2d) 
Jo) 
w ( A . 2d ) 

u(0) 
u ( R , 2d ) 

_ ( ° ) 
°(R.2d) 

5^ 

u ( R , 2d ) 

0.0 

0.424 

0.264 

0.623 

0.318 

0.256 

0.803 

0.050 

0.425 

0.264 

0.620 

0.334 

0.252 

0.755 

0.1C0 

0.428 

0.263 

0.614 

0.350 

0.250 

0.713 

0.150 

0.433 

0.262 

0.604 

0.366 

0.248 

0.677 

0.200 

0.439 

0.260 

0.593 

0.382 

0.247 

0.645 

0.250 

0.446 

0.259 

0.581 

0.398 

0.246 

0.618 

0.300 

0.454 

0.258 

0.569 

0.414 

0.246 

0.594 

0.350 

0.463 

0.258 

0.556 

0.430 

0.246 

0.574 

0.400 

0.473 

0.258 

0.545 

0.446 

0.248 

0.556 

0.450 

0.484 

0.258 

0.534 

0.462 

0.250 

0.541 

0.500 

0.495 

0.260 

0.524 

0.477 

0.252 

0.529 

0.550 

0.507 

0.261 

0.515 

0.493 

0.256 

0.518 

0.600 

0.520 

0.264 

0.508 

0.509 

0.259 

0.509 

0.650 

0.533 

0.267 

0.501 

0.525 

0.264 

0.502 

0.700 

0.547 

0.271 

0.496 

0.541 

0.269 

0.497 

0.707 

0.549 

0.272 

0.495 

0.543 

0.269 

0.496 

0.750 

0.561 

0.276 

0.492 

0.557 

0.274 

0.492 

0.800 

0.575 

0.281 

0.489 

0.573 

0.280 

0.489 

0.850 

0.590 

0.287 

0.486 

0.589 

0.286 

0.486 

0.900 

0.605 

0.293 

0.485 

0.605 

0.293 

0.485 

0.950 

0.621 

0.300 

0.484 

0.621 

0.300 

0.484 

1.000 

0.637 

0.308 

0.483 

0.637 

0.308 

0.483 
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the  standard  deviations  to  the  corresponding  averages,  varying 
from  about  0.5  to  0.7. 

To  correspond  to  customary  electrical  engineering 
practice,  we  now  examine  the  probability  law  of  1^°^  when 
represented  as  the  (cumulative)  probability  distribution  of 
its  logarithmic  value  in  decibels:  = 20  log^I^°*.  To 

represent  the  **esults,  we  use  a log-normal  graph  (a  log-normal 
distribution  [14]  is  a positively-sloped  straight  line  in  such 
a graph)  to  plot  the  probability  distributions  F{I^'  <_  i} 

[15] . These  plots  for  the  a = 0 and  a = 1 cases  are  snown  in 
Figure  4.  These  plots,  specified  by  (25)  and  (28),  are  not 
log-normal  and  therefore  are  not  straight  lines.  But  they 
are  approximately  linear,  up  to  "W5  percent  in  the  ( R , 2d)  and 
^60  percent  in  the  (A, 2d)  cases  for  a = 0 and  up  to  ^50  percent 
in  both  cases  for  a = 1;  and  they  flatten  to  one  at  higher 
percentiles.  The  flattening  at  the  higher  percentiles  is 
caused  by  the  induced  current  being  bounded  by  the  maximum 
response  of  the  basic  element  of  interaction  (when  mutual 
coupling  is  ignored) , which  provides  the  normalization  of  the 
curves.  The  effect  of  the  mutual  interactions  tends  to 
straighten  the  curves,  as  was  explained  by  (16)  and  its  fol- 
lowing reasoning,  and  is  shown  in  Figure  5.  Figure  5 plots 
the  mutual  coupling  included  distributions  F{I  = I ^ • |l  + G| 

< i},  with  mutual  coupling  strengths  o„  = 1 and  a_  = 30,  for 
the  corresponding  cases  in  Figure  4.  Excluding  the  extreme 
percentiles,  curves  in  Figure  5 clearly  exhibit  rather  close 
resemblances  to  straight  lines.  The  resemblance  is  better 
for  larger  a^. 

Within  this  resemblance,  I is  approximately  log-normal 
and  its  average  u and  standard  deviation  a are  related 
approximately  to  those  of  V = 20  log^I  by 
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a-1 , FOR  (A, 2d)  and  (R,2d)  CASES 

a =0,  FOR  (A, 2d)  CASE 

a=0,  FOR  (R,2d)  CASE 
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Figure  4.  Probability  Distribution  of  the  Normalized  Induced  Loop 

Currents  1^  for  the  Two-Dimensional  Cases  of  Equally 
Probable  Area  (A, 2d)  and  Equally  Probable  Radius  (R,2d), 
with  Size  Parameter  a =0,  and  a =1 
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Figure  5(a).  Probability  Distributions  of  the  Induced  Loop  Currents  I 
for  the  Two-Dimensional  Case  of  Equally  Probable  Area 
(A, 2d)  and  Equally  Probable  Radius  (R,2d)  with  Size 
Parameter  a =1  and  Mutual  Coupling  Strength  o„=l  and  ar= 30 
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Figure  5(b).  Probability  Distributions  of  the  Induced  Loop  Currents  I 
for  the  Two  Dimensional  Case  of  Equally  Probable  Area 
(A, 2d)  with  Size  Parametera=0  and  Mutual  Coupling 
Strengths  0=1  and  o =30 


INDUCED  LOOP  CURRENT 


Figure  5(c).  Probability  Distributions  of  the  Induced  Loop  Currents  I 
for  the  Two-Dimensional  Case  of  Equally  Probable  Radius 
(R,2d)  with  Size  Parameter  a=0  and  Mutual  Coupling 
Strengths  a =1  and  a =30 


= 30 


(31) 


20  I , 1 , T , . / a *2 ] | 

UY  " TrTTo-  |ln  U - j ln  [1  + (y)  J j 

°v  * ITTTo  Vln  t1  + O ] • 

Using  the  relation  (18)  that  links  the  o/p  of  I to  the 
o(o)/u(o)  of  I , we  immediately  have 


where 


20 

ln  10 


20 

ln  10 


(33b) 


(33c) 


The  superscripts  o and  00  of  Y represent,  respectively,  the 
cases  of  oG  + 0 and  o G -*•  °°.  The  two  limiting  standard  devia- 
tions of  the  logarithmic  currents, a - > and  a , > , are  plotted 

Y°°' 

in  Figure  6 as  functions  of  a^°^/u^°^.  The  shaded  area, 
representing  the  range  of  given  in  Table  2,  exhibits 

the  narrow  range  of  possible  values  for  the  corresponding 
log-normally  distributed  standard  deviations  a , , and  o , . . 

Y * y 00 

Figure  6 shows  a range  of  standard  deviations  for  Y encom- 
passing v4  to  ^8  dB  for  any  distribution-size  parameter  a, 
and  centered  around  ^6  to  v7  dB. 
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(0)*°  (°o)  DECIBLES 


and  oo  Mutual  Coupling  Strengths,  of  the  Logrithmic 


Currents  as  Functions  of  o^°V 
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Very  little  statistical  data  on  EM  coupling  to  large 
systems  is  available,  most  that  exists  was  taken  from  the 
PREMPT  program  [16].  These  tests  measured  the  maximum  peak- 
to-peak  current  response  to  broadband  pulsed  fields  and 
therefore  are  the  result  of  a more  complex  interaction  than 
the  CW  case  analyzed  in  this  paper;  however,  it  is  still 
interesting  to  see  if  the  analysis  and  the  data  follow  the 
same  general  trends. 

The  distributions  of  three  sets  of  such  measured  data, 
normalized  by  each  of  their  own  median  value,  are  shown  in 
Figure  7 [17].  They  are  rather  similar  to  each  other  and 
are  approximately  log-normal  with  standard  deviations  to 
their  logarithmic  values  at  ^5  to  ^7  dB.  A comparison  to 
the  theory  and  the  data  is  also  shown  in  Figure  7 by  overlaying 
the  median-value-normalized  theoretical  results,  for  the 
cases  of  (A, 2d)  and  (R,2d)  with  parameters  a and  a being 
0 and  1,  on  top  of  the  data.  The  comparison  shows  rather 
good  agreement  .for  the  ^40  percent  and  higher  percentiles, 
but  displays  a discrepancy  at  low  percentiles:  there  are 

fewer  low  (compared  with  the  median  value)  currents  in  the 
data  than  in  the  theory.  This  discrepancy  is  not  surprising. 

We  suspect  that  the  smaller  number  of  low  current  measurements 
is  caused  at  least  in  part  by  a systematic  bias  in  the  current 
data  that  would  result  from  an  engineering  judgment  to  avoid 
measuring  currents  that  are  obviously  uninteresting  because 
they  are  so  small  and  are  difficult  to  measure  in  the  background 
of  system  noise. 

To  the  extent  that  the  theoretical  curves  can  be  approxi- 
mated by  straight  lines  near  the  median  value  of  the  distri- 
butions, and  therefore  by  log-normal  distributions  over  this 
region,  the  standard  deviations  implied  by  their  slopes  are 
similar  to  those  of  the  data  (see  Figure  6)  . 


35 


II 

a 


ii 

a 


TD 

CM 


T3 

CM 


O 

II 


o 

II 

a 


o 
-a  ii 

g a 

TO 


O 

11 


T3 

C 

TO 


■o 

CM 


TD  ~0  ^ 
CM  CM  ^ 

< < 


*o 

CM 


"O 

CM 


iN^ano  Nvia3w/iN3ddno  w3d-oi->iV3d 


36 


PERCENTAGE  LESS  THAN  ORDINATE 
(NORMALIZED  2 SITE  INTERBAY  CURRENTS) 

Figure  7.  Peak-to-Peak  Induced  Current  Distribution  Data  for 
EMP  Simulations  Obtained  in  the  Program  PREMPT 
and  Overlay  of  2-D  Results 


3.2  THREE-DIMENSIONAL  ORIENTATIONS 


For  loops  oriented  equally  probable  in  the  three- 
dimensional  space,  as  might  roughly  be  the  case  for  a large 
electrical  facility,  we  can  obtain  the  induced  current 
probability  densities  by  the  same  procedures  used  in 
Section  3.1.  For  the  equally  probable  area  case, 


(loop) 

I (o)  (A, 3d) 
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0 < i < a 


a < i < 1 (34) 


and 
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The  three-dimensional  probability  densities  are  depicted 
in  Figure  8.  They  are  monotonical ly  decreasing  functions  of 
i,  as  expected  in  view  of  the  remark  made  at  (12)  and  (13) . 
Also,  they  favor  lower  current  values  more  than  their  two- 
dimensional  counterparts  do,  as  can  be  seen  easily  in  Figure  3 
or  by  observing  the  inequalities  in  (35)  to  (39).  Such  a 
behavior  is  intuitively  obvious  because  in  three-dimensions, 
unlike  in  two-dimensions,  there  is  more  solid  angle  within  a 
A9  near  the  perpendicular  to  a given  direction  than  there  is 
near  the  parallel.  Furthermore,  similar  to  the  two-dimensional 

cases,  we  have  u|R|3d)  1 3d)  and  °(r!  3d)  - ^ ( A ! 3d)  ' 


Table  3 lists  the  numerical  dependences  on  a of  the  averages 
and  standard  deviations.  They  exhibit  even  more  insensitivity 
to  the  changes  in  a than  the  two-dimensional  cases  do. 

Figure  9,  similar  to  Figure  4,  plots  in  log-normal  graphs 
the  distribution  of  for  the  three-dimensional  cases 

(A, 3d)  and  (R,3d)  with  a = 0 and  a = 1.  Finally,  we  point 
out  that  an  inclusion  of  the  mutual  coupling  effect,  similar 
to  that  in  Section  3.1,  again  renders  a cumulative  current 
distribution  that  is  in  its  central  part  roughly  log-normal 
and  results  in  standard  deviations  for  20  log^Q  1 even  slightly 
more  closely  packed  near  ^ 6 dB . The  last  statement  can  be 
easily  seen  by  reading  the  a^°^/u^°^  values  of  Table  3 into 
Figure  6.  A comparison  of  the  PREMPT  data  and  the  three- 
dimensional  loop  distribution  results,  similar  to  the  two- 
dimensional  results  shown  in  Figure  7,  is  shown  in  Figure  10. 

As  before,  the  trends  in  the  theory  and  the  data  correspond 
quite  well. 
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Figure  8.  Probability  Densities  of  Induced  Currents  in  Loops  Uniformly 

Oriented  in  3-Dimensional  Space,  P^m^  (1)  and 
Hoop)  I u (A ,3d) 

P (0),  (i ),  Respectively,  of  Equation  34  fbr  Equally  Probabl 
I ( R , 3d ) 

Loops  Areas  and  of  Equation  37  for  Equally  Probable  Loop's 


TABLE  3.  AVERAGES  AND  STANDARD  DEVIATIONS  FOR  MAGNETIC  LOOPS  EQUALLY 
PROBABLY  ORIENTED  IN  THREE-DIMENSIONAL  SPACE 
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NORMALIZED  INDUCED  LOOP  CURRENT 


Figure 9 Probability  Distribution  of  the  Normalized  Induced  Loop 

Currents  1^  for  the  Three-Dimensional  Cases  of  Equally 
Probable  Area  (A, 3d)  and  Equally  Probable  Radius  (R,3d) 
with  Size  Parameter  a =0 , at =1 
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Figure  10.  Peak-to-Peak 
EMP  Simulati 
and  Overlay 


3 . 3 REMARKS 


The  results  of  the  model,  taken  together  with  the 
limited  available  data,  lead  to  the  following  tentative 
conclusions.  First,  our  statistical  coupling  models  are  not 
of  the  type  that  have  a large  number  of  parameters  and  can  be 
made  to  fit  virtually  any  data  by  adjusting  these  parameters. 

On  the  contrary,  our  models  have  very  few  parameters  and 
yield  distributions  for  their  major  central  part  insensitive 
to  the  values  of  those  few  parameters  (see  Section  3.2  for 
three-dimensional  cases).  In  view  of  this,  the  correspondence 
between  the  data  and  the  models,  as  pointed  out  in  the  previous 
sections,  may  have  important  implications  in  two  aspects.  On 
the  one  hand,  the  EM  coupling  to  large  systems,  when  dominated 
by  low  frequency  magnetic  fields,  is  largely  insensitive  to 
the  coupling  detail  and  yields  a distribution  whose  central 
part  is  nearly  log-normal  with  a standard  deviation  of  about 
6 dB.  On  the  other  hand,  precisely  because  of  the  insensi- 
tivity of  the  major  central  part  of  the  distribution  to  the 
model,  one  may  not  be  able  to  use  it  to  determine  the  detailed 
parameters  of  the  coupling  (e.g. , the  loops'  orientations, 
their  relative  sizes,  their  mutual  coupling  strengths,  etc.). 
Second,  the  shapes  of  the  extreme  percentiles  of  the  distribu- 
tion may  depend  on  and  be  sensitive  to  the  detailed  nature 
of  the  coupling.  This,  combined  with  the  first  conclusion, 
would  make  apparently  plausible  assumptions  concerning  the 
data  distribution  function  (e.g.,  log-normal)  and  the  extrapo- 
lation to  the  extreme  percentiles  based  on  the  central  part 
of  the  distribution  susceptible  to  substantial  errors.  Con- 
sequently, any  conclusions  dependent  upon  such  assumptions 
and  extrapolations  might  be  susceptible  to  substantial  errors. 
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where  the  last  inequality  uses  the  assumption  that  the 
thickness  t of  the  loop  wire  is  much  larger  than  the 
conducting  loops'  skin  depth.  Using  this,  (1)  implies 
(2)  . 

8.  See  any  textbook  on  probability  theory,  e.g.,  E.  Parzen, 
Modern  Probability  Theory  and  Its  Applications,  Chapter  7 
(Wiley  & Sons,  i960) . 

9.  For  thin  short  wires,  linear  wire  current  distribution 
gives  f ( M , L ) = M/(L/2),  and  D <<  L gives  the  static  self- 
capacitance 
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Again,  the  last  inequality  uses  the  assumption  that  the 
wires'  diameters,  although  small  compared  with  the  wires' 
lengths,  are  not  small  compared  with  the  conductor's 
skin  depth  6W  at  the  frequencies  considered.  Using  these, 
we  obtained  (6)  and  (7)  . 

10.  W.  Feller,  An  Introduction  to  Probability  Theory  and  Its 
Appl ication , Vol.  2,  pp.  32  (Wiley  & Sons,  1971). 

11.  This  can  also  be  shown  formally  by  using  circuit  analysis 
notations.  For  a linear  network,  we  have  for  the  kth 
circuit 
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See  any  probability  textbook,  e.g.,  Reference  8,  Vol.  1, 
Chapter  10  (third  edition,  1967).  Strictly  speaking,  the 
fields  Ej  from  all  other  dipoles  j / i are  not  totally 
independent  of  each  other,  because  of  their  mutual  couplings. 
Intuitively,  these  Ej's  are  approximately  and  essentially 
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independent  random  variables  since  their  deterministic 
mutual  coupling  mechanisms  are  controlled  by  the  independent 
random  variables  associated  with  their  geometric 
configuration.  However,  we  have  not  been  able  to  prove 
such  an  approximate  and  intuitive  independence  and  quantify 
rigorously  its  implication  to  the  applicability  of  the  central 
limit  theory. 

13.  Use  the  integration  results 
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14.  See,  e.g.,  J.  Aitchson  and  J.  A.  Brown,  The  Log-Normal 
Distribution , The  Cambridge  University  Press  (Reprinted 
1963,  first  printed  1957). 
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ing densities  as: 
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APPENDIX  A 


SPECIAL  CASES  FOR  ELECTRIC  WIRES 

The  results  of  the  short  electric  dipole  coupling 

probability  distributions  of  I ^ for  several  interesting 

special  cases  are  given  in  this  section.  Here,  the  nota- 

2 2 

tions  used  are:  K = £ In  b/(b  In  a)  <_  1;  x^(i)  the  solution 

2 2 

of  x^  In  b = ib  In  x^;  and  1 <<  a _<  L/D  = X _<  b.  The  analysis 
follows  procedures  similar  to  those  in  Section  3,  but  with 
the  basic  coupling  described  by  (6)  to  (10). 

A . 1 COPLANAR  ORIENTATIONS 

For  the  short  wires  all  lying  in  the  9 = 7r/2  plane  (see 
Figure  2)  and  oriented  equally  probable  in  all  azimuthal 
directions,  a plane  wave  incident  from  the  +z  axis  induces 
the  following  current  probabilities. 

When  both  wires'  lengths  X and  the  relative  position  of 
sampled  current  locations  S are  fixed,  we  obtain  for  the 
normalized  current  1^°' 
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Here  the  subscripts  (l,2d)  indicate  the  first  case  considered 
in  two-dimensions.  It  has  a uj°^2d)  = = and  a 

0 ( l!  2d)  = (1/2  " V*2)172  = 0-31. 

When  X is  equally  probable  in  the  interval  [a,b]  and  S 
is  fixed,  we  obtain 
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Its  average  and  standard  deviation  can  be  obtained  numerically. 

When  X is  fixed  and  the  current  sampling  loc^ions  are 
equally  probable  anywhere  along  the  wire,  S uniform  in  [0,1], 
we  obtain 
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which  has  a u ^ ° ^ 2d ) = ^‘^TT  = and  a °(3^2d)  = " l/^2)^2 

= 0.26. 

To  investigate  the  effects  of  nonnormal  incidences,  we 
let  the  incidence  make  an  angle  0q  i-  0 with  the  -z  axis. 

Now  the  polarization  angle  f influences  the  results.  For 
example,  if  f is  a fixed  angle  </>  and  both  X and  S are  fixed, 
then 
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where  9'  is  the  angle  between  the  z axis  and  the  incident 
o ~ 

polarization  E,  cos  9'  = cos  b sin  9 . The  (A-4)  has  a 
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where 


f (;)  = 4 
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Motice,  of  course,  both  (A-5)  and  (A-6)  reduce  to  (A-l)  as 
9q  - 0 . The  nature  of  the  probability  densities  (A-l)  to  (A-5) 
is  depicted  by  Figure  A-l. 

A. 2 THREE-DIMENSIONAL  ORIENTATIONS 

For  wires  with  equally  probable  orientations  in  a three- 
dimensional  space,  the  induced  currents’  probability  distri- 
butions are  independent  of  the  incident  and  the  polarization 
directions.  Except  the  trivial  case  of  fixed  X and  fixed  S, 
which  incurs  a uniform  currents'  probability  distribution, 
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Figure  A-l. 


Probability  Densities 
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i ) of  Induced  Currents 


on  Short  Electric  Wires  Uniformly  Oriented  in  a Plane 


several  interesting  cases  are  given  in  the  following  and  their 
resulting  probability  densities  are  sketched  in  Figure  A-2. 

If  X is  uniform  in  [a,b]  and  S is  fixed,  then 


(wire) 

PI(o) (1,3d) 


(i) 


b^ln  x 

2.  . 
x In  b 


dx 


b^ln  x 

2 , . ' 
x In  b 


0 < i < K 


K < i < 1. 

(A— 7) 


If  X is  fixed  and  S is  uniform  in  [0,1] , then 
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which  has  an  average  1/4  and  a standard  deviation  /1/9  - 1/16 
= 0.22.  If  both  X and  S are  uniformly  distributed,  respectively, 
in  [a,b]  and  [0,1],  then 


, -b  .1 

^ , i f In  x j , / di’ 

— dX+/  T- 


(wire) 

PI(0)  (3,3d) 


(i)  = 


(b  - a)  In  b 


/ 


In  x 


dx , 0 < i < k 


x1(i’) 


r 


In  x 


i x : ( i 1 ) 


dx , K < i < 1 


(A-9) 


A- 5 


(wire) 


Figure  A-2.  Probability  Densities  of  Induced  Currents 


on 


Short  Electric  Wires  Uniformly  Oriented  in  3-Dimensional 
Space,  Except  for  Eq.  ( A- 1 0 ) which  is  for  all  Angular 
Orientations  Being  Fixed. 


A-6 


Now,  to  distinguish  the  separate  effects  as  caused  by  the 
randomness  of  the  lengths  X and  that  of  the  relative  positions 
S,  we  consider  the  wires'  orientation,  as  well  as  the  incident 
wave  direction  and  polarization,  to  be  fixed.  Then,  a uniform 
X in  [a,b]  and  a fixed  S gives 


(wire) 

PI(o) (4,3d) 


(i)  = 


[bln  xl  (i) 


(b  - 


"a)  x^  (i ) 


In  b ( 2 1 n x ^ ( i ) 


K < i < 1. 


- i; 


(A- 10 ) 


However,  a fixed  X and  a uniform  S in  [0,1]  gives  simply  a 
uniform  I ^ in  [0,1],  as  it  should  be.  Finally,  the  proba- 
bility density  for  both  the  X uniform  in  [a,b]  and  the  S 
uniform  in  [0,1]  is  exactly  (A-7)  again.  From  these,  we  clearly 
see  that  as  far  as  the  induced  currents'  probability  is  con- 
cerned, the  effect  of  randomness  in  the  current  measurement 
position  S,  uniform  or  fixed,  is  the  same  as  that  of  randomness 
in  the  wire  orientations  in  three-dimensional  space. 


A-7 


\ 


APPENDIX  B 


ELLIPTIC ALLY  POLARIZED  INCIDENCE 


If  the  incident  monochromatic  EM  wave  is  elliptically 
polarized,  then  the  coherence  of  the  two  superposing  perpen- 
dicular  linear  field  components  substantially  changes  the  nature 
of  the  probability  density  of  the  basic  induced  current  1^°^. 

To  bring  out  this  effect,  we  shall  consider  a simple  case  of 
an  elliptically  polarized  wave  incident  upon  a short  wire  of 
direction  L from  +z  axis,  with  a semi-major  axis  of  polarization 
in  the  y-direction  and  an  angle  of  ellipticity  x in  [ 0 , tt/4  ] 
inc  inc 

defined  by  |E  /E^  i tan  x £ 1 • Notice  that  by  definition 

the  x 0 case  is  a linear  polarization,  and  the  x "*■  one 

is  a circular  polarization  [19].  It  follows  then  that 


With  this  £ replacing  the  n in  (3)  and  (6),  we  obtain  the 
following  results. 

B.l  COPLANAR  ORIENTATIONS 

When  the  wires  lie  in  the  9 = tt/2  plane  but  orient 
themselves  with  equal  probability  azimuthally  and  have  fixed 
X and  S,  we  obtain 

0 £ i £ tan  x 

tan  x £ i £ 1 

( B—  2 ) 

B-l 


This  reduces  to  (A-l)  when  x + 0 for  a linear  polarization,  as 
it  should,  becomes  (i  - 1)  when  x ■*  ir/4  for  a circular  polar- 
ization which  does  not  distinguish  one  direction  in  the  plane 
from  another,  and  has 


u 


(o) 

X ( 1 , 2d ) 


e(  J,  V 1 - tan2x  ) 


X - 0 

► 


2/tt 


X -*•  ^ ^ ( B — 3 ) 


a {o) 
Xd»2d) 


sec 


u 


(o) 

X ( 1 » 2d ) 


x - 0 1 

\2 


X - tt/4  0 ( B—  4 ) 


where  E(<$,k)  is  the  elliptical  integral  of  the  second  kind. 

Table  B-l  shows  and  a'!0'-,,  as  functions  of  .. 

(X»2d>  ( x » 2d) 


TABLE  B-l.  AVERAGES  AND  STANDARD  DEVIATIONS  FOR  INDUCED  CURRENT  ON  WIRES 
UNIFORMLY  ORIENTED  IN  A PLANE,  WITH  AN  ELLIPTICALLY  POLARIZED 
AND  NORMALLY  INCIDENT  PLANE  WAVE 


X 


(degree) 

tan  x 

u 

0 

o/ u 

0 

0.0 

0.637 

0.308 

0.483 

4.5 

0.079 

0.643 

0.299 

0.466 

9.0 

0.158 

0.659 

0.280 

0.425 

13.5 

0.240 

0.679 

0.260 

0.382 

18.0 

0.325 

0.707 

0.231 

0.327 

22.5 

0.414 

0.738 

0.203 

0.276 

27.0 

0.510 

0.773 

0.178 

0.230 

31.5 

0.613 

0.818 

0.136 

0.166 

36.0 

0.727 

0.868 

0.100 

0.115 

40.5 

0.854 

0.928 

0.066 

0.071 

42.5 

0.916 

0.959 

0.025 

0.026 

45 

1.000 

1.000 

0 

0 
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B . 2 THREE-DIMENSIONAL  ORIENTATIONS 


When  the  wires  orient  equally  probable  in  the  three 
dimensional  space,  the  probability  density  for  a uniform 
in  [a,b]  and  a fixed  S is 


X (1.3d) 


(i)  -V 


b u^(xi 

/ lx  J du 

a u (x) 

o 

« i bur 

/ %»  / y1 


U (X) 


/ d.  / 

» '4  f X I 

0 

1 

y"  dx  y*  du  * y*  ix  y-  du 

4 fx) 


u,  :*) 

r 

0 ( X ) 


Z1*  / 


(.  I 


, 9 « l < K tan 


. K tan  » i i mn  IK,  tan  <) 
2 


2iu  1-4^  \ 
‘ x In  o / 


1 


» ( a)  ^l-u^[u4-j^x)  I x)  *u  ^ tan  4 ( 


. tan  t « 1 < K 


K I 1 i tan  j 


, sax  ; K . tan  r ) « i.  • l 


which,  of  course,  reduces  to  (A-7)  as  x + 0.  Here,  u (x) 

2 2 2 2 
ib  In  x/(x  In  b) ; u^(x)  = ib  In  x/(tan  x * x In  b) ; uq(x^ 

from  which  x.  - b as  i -*•  1 and  x,  -►  a as  i -*•  K;  and  u,  (x 
1 1 1 c 

from  which  x -*•  b as  i - tan  x and  x -*■  a as  i -*•  K tan  x* 
o A o 

Several  special  cases  of  (B-5)  are  worth  noticing. 

wires  with  the  same  lengths,  a -*■  b and  (B-5)  reduces  to 


(B-5) 

) i 1 
) = 1 

For 


B-3 


(wire , x) 


, 0 < i < tan  X 


./(I-i2)/(l-tan2x) 


dv • G ( v ; i , X ) 


, tanx  < i < 1 


where 


G ( v ; i , X ) 2 J 


^1-v2) [1-v2 (l-tan2x) ] [l-i2-v2 (l-tan2x) ] 


This  density  further  reduces  to  a uniform  one  as  x -0  and  to 

i//l-i2 , with  an  average  n/4  v 0.79  and  a o = (2/3-tt2/16  ) 1,/2 
^ 0.22,  as  x tt/4  . In  general,  the  (B-6)  starts  with  value 
0 at  i = 0 and  ends  with  value  ( l-tan2x)  1//2  at  i = 1. 

The  general  behavior  of  ( B—  5 ) starts  with  0 at  i = 0, 
reaches  a maximum  in  0 < i < 1,  and  ends  with  0 at  i = 1. 
Further,  as  x ■*  it/4,  (B-5)  reduces  to 


xL(i) 


i 

b - a 


b In  x 
x2ln  b 


0 < i < K 


-o'-) 


K < i < 1 


B-4 


From  the  above,  we  see  that  for  spatially  randomly 
oriented  wires  with  no  preferred  direction,  the  elliptically 
polarized  incidence,  as  compared  to  a linearly  polarized  one, 
makes  it  less  probable  for  smaller  (near  zero)  induced 
currents  and  pushes  the  most  probable  current  to  a higher 
value  (toward  the  maximum  possible  current) . A basic  differ- 
ence from  the  linear  polarization  case  is  that  the  probability 
densities  are  not  a monotonically  decreasing  function  of  the 
induced  current  amplitude.  In  view  of  the  similar  nature 
with  which  the  polarizations  enter  into  the  basic  coupling 
mechanisms,  (1)  and  (6),  these  conclusions  apply  to  the 
magnetic  coupling  cases  as  well  as  to  the  electric  ones. 
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